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1 Introduction 



1.1 Formulation of the problem 

The nonlinear second-order evolution equation 

ut = {f{u)u^)^ + {g{u)uy)y + {h{u)u:,)^ (1.1) 

describes the heat conduction in anisotropic materials whose physical char- 
acteristics such as the thermal conductivity are affected by the temperature. 
Physical applications and interesting mathematical properties of Eq. f ll.ip and 
of its extension to the case of existence of an external source q{u), 

Ut = {f{u)ux)x + {g{u)uy)y + {h{u)u^)z + q{u) (1.2) 

are discussed in pD| (see also [2], Section 10.9). 

The functions f{u),g{u), h{u) are positive according to their physical mean- 
ing. So, we consider Eqs. f ll.ip and 01.21] with arbitrary positive coefficients 
f{u), g{u), h{u). Furthermore, we will assume that these coefficients are lin- 
early independent and that none of them is constant, i.e. 

f'iu)^0, g'iu)^0, h'iu)^0. (1.3) 

Note, that Eq. (11.11) has a conservation form whereas Eq. (11.21) with 
q{u) does not have such form. A conservation form is useful in many 
respects, e.g. in qualitative and numerical analysis. Moreover, possibility of 
different conservation forms can be helpful. Therefore we will construct various 
conservation laws for Eq. (II. ip using the method of nonlinear self-adjointness 
[3] and investigate the question on existence of conservation laws for Eq. (II. 2p 
with specific values of the source term q{u) ^ 0. 

1.2 Nonlinear self-adjointness 

Recall the definition of nonlinear self-adjointness. Let us consider a second- 
order partial differential equation 

F(a;,'u,'U(i),'U(2)) = 0, (1.4) 

where u is the dependent variable, and M(2) are the sets of the first-order 
partial derivatives Ui and the second-order derivatives Uij of u with respect to 
the independent variables x = [x^, . . . , x"). The adjoint equation to Eq. (11.40 
is 

F*{x,U,V,U(l),V(l),U(^2),V(^2)) =0, (1.5) 
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where F* is defined by 

F [x,U,V,U(i),V(i),U(^2),V(2)) = ■ (1.6) 

Here f is a new dependent variable and f(i),f(2) are the sets of its partial 
derivatives. Furthermore, 6{vF)/6u denotes the variational derivative of vF : 

6{vF) _ d{vF) (divFY^ ^ ^^^^ fd{vF) 



du ou V oui J V ouik 

where the total differentiations are extended to the new dependent variable v : 

^ d d d d d ,^ ^, 

Di = — + Ui— + Vi— + Uij- VVij- \ . (1.7) 

OT* OU OV OUj OVj 

Eq. (11. 4p is said to be nonlinearly self-adjomt [3J if the adjoint equation 
(II. 5p is satisfied for all solutions u of the original equation (ll.4p upon a sub- 
stitutional 

v = ip{x,u), f 0. (1.8) 

The condition that the function ip that does not vanish is significant. The 
condition for the nonlinear self-adjointness can be written in the form 

F*{x,U,(p,U(^l),ip^l),U(^2),^{2)) = AF(x,M,M(i),'U(2)), (1.9) 

where A = A(x, u, . . .) is an undetermined variable coefficient and ip(^2) 
denote the derivatives of the function (p. E.g. is the set of the first-order 
total derivatives 



Eq. ( 11. 9p should be satisfied identically in all variables 



(2)- 



1.3 Conserved vector associated with symmetries 

The general result on construction of conserved vectors associated with sym- 
metries of nonlinearly self-adjoint equations demonstrated in [3] leads to the 
following statement for the second-order equation (11.40 . 

Let (11.40 be nonlinearly self-adjoint and admit a one-parameter point trans- 
formation group with the generator 

^In general, the substitution (jl.8l) can be of the form v = (^(x, m, U(-x)) . 
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Then the vector 



C = W 



dui \ du. 



is a conserved vector for Eq. (11.41) . i.e. satisfies the conservation equation 
Here 

W = ri- i^Uj (1.13) 
and C is the formal Lagrangian for Eq. (I1.4p given by 

C = vF. (1.14) 

It is assumed that the variable v and its derivatives are ehminated from the 
right-hand side of Eq. (11.111) by using the substitution (II. 8p . where the function 
ip{x,u) is found by solving the nonlinear self-adjointness condition (II. 9p . 

2 Investigation of nonlinear self-adjointness 
2.1 Substitution ( ITTSi ) for Equation ( OTB 



Since Eq. (11.11) has the conservation form (11.120 . it is nonlinearly self-adjoint 
by Theorem 8.1 from p|. Let us find the corresponding substitution (II. Sp . 
We write Eq. ( 11. ID in the form (\1A\\ : 

F = -ut + f{u)u,,+g{u)uyy + h{u)u,, + f'{u)ul + g'{u)ul + h'{u)ul = 0, (2.1) 

insert the expression for F in (11.61) and after simple calculations obtain the 
following adjoint equation (II. 5p to Eq. (II. ip : 

F* = vt + f{u)v^^ + g{u)vyy + h{u)vzz = 0. (2.2) 

In our case the substitution (ll.Sp has the form 

V = ip{t,x,y,z,u). (2.3) 

Its derivatives are written 

Vt = Dt{(f) = (fuUt + ft, = D^{(f) = i^uUx + V^x, 



V. 



y 



Dy{lf)=lfuUy + !fy, V ^ = D ^ {(f) = If ^ + (fi ^ 



V^x = Dl{(f) = IfuUxx + (fuuUl + 2(fxuUx + fxx, (2.4) 

= ^liv) = 'fuUyy + fuuUl + IfyuUy + fyy, 

V^^ = Dl{if) = ifuUzz + VuuUl + 2(pzuUz + (fzz. 
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Now we take the nonlinear self-adjointness condition fll.9p . where F and F* 
are given by fl2.ip and fl2.2p . respectively: 

vt + f{u)v^^ + g{u)vyy + h{u)vzz (2.5) 
=A [-Mt + f{u)uxx + g{u)uyy + h{u)uzz + f{u)ul + g'{u)ul + 

where the corresponding derivatives of f in the left-hand side should be re- 
placed with their expressions (12.40 . First we compare the coefficients for Ut in 
both sides of Eq. (12.51) and obtain 

A = -ipu- 

Then the coefficients for Uxx,Uyy,Uzz yield: 

f{u)ipu = -f{u)ipu, g{u)(fu = -g{u)ipu, h{u)ifu = -h{u)ifu- 

By our assumption, the functions f{u),g{u), h{u) do not vanish. Therefore the 
above equations yield that = 0. Hence, ip = ip{t,x,y, z) and therefore 

A = 0, Vt = (ft, Vxx = ^xx-i "Vyy = ^yyi '^zz = ^zz- 

Then Eq. (12. 5p becomes 

(ft + f{u)ipxx + g{u)<^yy + h{u)ipzz = 0. (2.6) 

Since f{u),g{u),h{u) are linearly independent and obey the conditions (II. 3p . 
whereas ip does not depend on u, Eq. (12.60 yields 

ft = 0, (fxx = 0, (Pyy = 0, iPzz = 0. (2.7) 

The general solution of Eqs. (12. 7p is given by 

ip = ai xyz + a2xy + a^xz + a^yz + a^x + aQy + ttY z + as 

with arbitrary constant coefficients ai, . . . og. This proves the following. 

Proposition 2.1. Eq. (11.10 satisfies the nonlinear self-adjointness condition 
(11.91) with the substitution (12.30 of the form 

V = ai xyz + a2xy + a^xz + a4yz + a^x + Gey + aj z + Og. (2.8) 
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2.2 Two-dimensional equation with a source 

Let us consider Eq. (11 ■2p , for the sake of simplicity, in the case of two spatial 
variables x,y : 

Ut = {f{u)ux)x + {g{u)Uy)y + q{u). (2.9) 
The adjoint equation has the form 

F* =vt + f{u)v.^^ + g{u)vyy + (^{u)v = 0. (2.10) 

Repeating the calculations of Section 12.11 we obtain the following equation 
for the nonlinear self-adjointness of Eq. (12. 9p (compare with Eq. (12. 6p ): 

ipt + f{u)ip^^ + g{u)ifyy + q'{u)if = 0. (2.11) 

If f{u),g{u) and q{u) are arbitrary functions, Eq. (12. lip yields (compare with 
Eqs. HMD): 

V^i = 0, (p^^ = 0, Lpyy = 0, Lp = 0. (2.12) 

These equations show that a substitution of the form (11.80 does not exist. 
Indeed, the last equation in (12.120 contradicts the condition f 0. Hence, Eq. 
(12.90 with the arbitrary source q{u) is not nonlinearly self-adjointness with the 
substitution of the form (11.80 . 

However, Eq. (12. 9p with sources of particular forms can be nonlinearly 
self-adjoint. For example, let 

q'{u) = rf{u), r = const. (2.13) 

Then Eq. (12.110 becomes 

ipt + fiu)[ip^^ + TLp] + g{u)iPyy = 

and yields (compare with Eqs. (12.120 ): 

<pt = 0, ipyy = 0, (f^^ + rip = 0. (2.14) 

The solution to Eqs. (12.140 has the form 

ip = a{x)y + b{x) , (2-15) 

where a{x) and b{x) arbitrary solutions of the linear second-order ODE 

w" + rw = 0. (2.16) 

Eq. (12.130 shows that the source strength increases together with the tem- 
perature, i.e. q'{u) > 0, if r = > 0, and decreases, q'{u) < 0, if r = —5^ < 0. 
Having this in mind and denoting 

jr(^) = j f{u)du (2.17) 
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we consider two particular forms of Eq. fl2.9p : 

ut = {f{u)u,j),j; + {g{u)uy)y + uj'^J^{u), uj = const., (2.18) 

and 

Ut = {f{u)u^)^ + {g{u)uy)y - 6'^I'{u), 5 = const. (2.19) 
In the case fl^T^ Eq. fIZTB]) is written 

w" + uP'w = 

and yields 

w = Ci cos(a;a;) + C2 sin(a;x). 

Hence 

a(x) = Ai cos{ujx) + A2 sin(ci;x), b{x) = Bi cos(a;a;) + B2 sin(ci;x) 

with arbitrary constants Ai, A2, Bi, B2. We we substitute these expressions in 
Eq. fl2.15p and arrive at the following statement. 

Proposition 2.2. Eq. (12.181) satisfies the nonlinear self-adjointness condition 
(11.91) with the substitution (12.31) of the form 

V = {Aiy + Bi) cos(u;x) + (A2I/ + B2) sin(a;x). (2.20) 

In the case (l2A9l) Eq. (IXTHD is written 

w" -5'^w = Q 

and yields 

w = Cie^^ + Cae--^^ 
Proceeding as above we arrive at the following statement. 

Proposition 2.3. Eq. (12.191) satisfies the nonlinear self-adjointness condition 
(11.91) with the substitution (12.31) of the form 

V = {Aiy + 5i) e-^^' + (A2I/ + B2) (2.21) 

2.3 Remark on materials with specific anisotropy 

The situation is different if the conditions (11.31) are not satisfied. Let, e.g. g{u) 
be a positive constant, g = k. Then Eq. (II. ip has the form 

Ut = (/(u)mx)x + kuyy + {h{u)u^)^. (2.22) 
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In this case Eqs. fl2.7p are replaced by the following equations: 

ipt + kifyy = 0, Lf^^ = 0, ip^^ = 0. (2.23) 
The second and third equations of the system fl2.23p yield 

ip = a{t, y)xz + l3{t, y)x + 7(t, y)z + a{t, y). 

The first equation f l2.23p shows that a{t, y), I3{t, y),'y(t, y) and a(t, y) solve the 
adjoint equation 

vt + kvyy = (2.24) 

to the linear heat equation 

Ut — kUyy = 0. (2.25) 
Thus, we have demonstrated the following statement. 

Proposition 2.4. Eq. fl2.22p satisfies the nonlinear self-adjointness condition 
(11. 9P with the substitution (12. 3p of the form 

V = a{t, y)xz + Pit, y)x + -f{t, y)z + a{t, y), (2.26) 

where a(t,y), /3(t,y),'j{t,y) and <j(t,y) are any solutions of the adjoint equa- 
tion (I2.24P to the linear heat equation (I2.25p . 

Combining Propositions 2.2 and 2.3 with Proposition 2.4, we obtain the 
following statements. 

Proposition 2.5. The equation 

Ut = {f{u)u^)^ + kuyy + u^I'i^u), f{u) = F'{u), (2.27) 

satisfies the nonlinear self-adjointness condition (11. 9p with the substitution 
()2.3p of the form 

V = a{t, y) cosiujx) + /3(t, y) sin(c<;x), (2.28) 

where a(t, y) and /3(t, y) are any solutions of the adjoint equation (I2.24p to the 
linear heat equation (12. 25 p . 

Proposition 2.6. The equation 

ut = {f{u)u.,),, + kuyy-5^j^{u), f{u) = 7'{u), (2.29) 

satisfies the nonlinear self-adjointness condition (II. 9p with the substitution 
flO) of the form 

V = a{t,y)e'^ + l3{t,y)e-'\ (2.30) 

where a(t, y) and /3(t, y) are any solutions of the adjoint equation (I2.24p to the 
linear heat equation (I2.25p . 
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3 Conservation laws 



3.1 Computation of conserved vectors for Equation (11.11) 

Here we construct the conserved vector (11 .lip for Eq. (11. ip . 



ut = f{u)uxx + g{u)uyy + h{u)u^^ + f{u)ul + g'{u)ul + h'{u)u 



(3.1) 



associated with its translational symmetries. We specify the notation by writ- 
ing the sjTiimetry generator (11.101) in the form 



d 



d 



dt dx 
Then the expression (I1.13P becomes 



dy 



dz 



+ r] 



d_ 

du 



(3.2) 



w = r^- eut - - euy - e 



and the conservation equation (I1.12p means that the following equation holds 
on the solutions of Eq. (13.11) : 

Dt{C')+Dx{C^)+Dy{C')+D,{C')=0. (3.4) 

The formal Lagrangian (11.14p for Eq. ( 13. ip is 

C = v[f{u)uxx + g{u)uyy + h{u)uzz + f'{u)ul + g'{u)ul + h'{u)ul - ut]. (3.5) 

Due to the specific dependence of the formal Lagrangian (13. 5p on the deriva- 
tives Ui,Uij, the components of the vector (II. lip are written 

dC 




Substituting here the explicit expression (13. 5p of C we obtain: 

Ci = -Wv, 

C^ = W [f\u)uxv - f{u)vx] + f{u)vDx{W), 
C'' = W {g'{u)uyv - g[u)vy\ + g{u)vDy{W), 
C^ = W [h\u)u,v - h{u)v,] + h{u)vD,{W). 



(3.6) 
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Eqs. fll.ip and fll.2p with arbitrary coefficients f{u),g{u), h{u) are invariant 
under the groups of translations of t, x, y, z with the generators 

d d d d 

-^1 = TTT , -^2 = 7^, -^3 = 5 X4 = — • (3.7) 

at ax ay az 

Eq. (11.11) has also a dilation symmetry. Moreover, both equations (II. ip and 
(11.21) may have more symmetries in certain particular cases [2], but we don't 
consider them here. 

Let us apply the formula (13. 6p to the symmetry X2. The corresponding 
quantity (13.30 equals W = —Ux- We substitute it in (13. 6p and obtain 



= vux, 

= -f'{u)vul + f{u)uxVx - f{u)vu^ 

= -g'{u)vUxUy + g{u)uxVy - g{u)vuxy, (3.8) 

= —h' {u)vUxUz + h{u)UxVz — h{u)vUxz- 

We have to substitute here the expression (12. 8p for t>, 

V = ai xyz + a2xy + a^xz + a^yz + a^x + a^y + z + ag. (12.80 

Since f is a given function whereas u is any solution of Eq. (13. ip . we want 
to simplify the conserved vector ( 13.80 by transforming it in an equivalent con- 
served vector which conserved density contains u instead of Ux- To this end, 
we use the identity vux = Dx{uv) — uvx and write in (13. 8p in the form 

=C^ + Dx{uv), 

where ^ 

= -uVx. (3.9) 

Then we transfer the term Dx{uv) from to using the usual procedure 
(see, e.g. [3], Section 8.1). Namely, since the total differentiations commute 
with each other, we have 

Dt{C' + Dx{uv)) + Dx{C^) =Dt{C')+Dx{C^ + Dt{uv)). 

Therefore the conservation equation (13.41) for the vector (13.80 can be equiva- 
lently rewritten in the form 

Dt{C')+Dx{C^)+Dy{C^)+Dz{C^)=0, (3.10) 

where has the form (13. 9p and is given by 

= C^ + Dt(uv). 



11 



Let us work out the above expression for C^. Invoking that Dt{v) = due to 
Eq. (gS]), and using Eqs. ([HI]), ([XI]) we have 

C^ = C^ + vut = C^ + v [f{u)u,, + f\u)ul + Dy {g{u)uy) + {h{u)u,)] . 

We substitute here the expression of from Eqs. (13 .Sj) and obtain: 

& = f{u)u^v^ + vDy[g{u)uy) +vD^[h{u)u^). 

We simphfy the latter expression for by noting that 

vDy[g{u)uy) = Dy[g{u)vuy) - g{u)uyVy, 

vDz[h(u)uz) = Dz[h(u)vUz) — h{u)uzVz 

Therefore we transfer the terms Dy[g{u)vUy) and Dz{h{u)vu^ to and C^, 
respectively, and obtain: 

= f{u)u^v^ - g{u)uyVy - h{u)uz. (3.11) 
Now the components and of the vector (13. 8|) become: 

= + D^{g{u)vuy), C^ = C^ + D,{h{u)vuz). 
After substituting here the expressions of and C'^ from (13. 8|) we have 

= g{u) {uxVy + UyV^) , = h{u) [u^v^ + u^v^) . 

Combining these expressions with (13.91) . (13.111) and ignoring the tilde, we arrive 
at the following conserved vector which is equivalent to (13.81) : 

= f{u)u^v^ - g{u)uyVy - h{u)UzVz, (3.12) 
= g{u){u.^Vy + UyV^), 

The vector (I3.12p involves the first-order derivatives of the variable v given by 
Eq. (12. 8p . Therefore the vector (13.121) contains seven parameters oi, . . . , 07. In 
fact, it is a linear combination of seven linearly independent conserved vectors 
obtained from (I3.12p setting by turns one of the parameters a, equal to 1 and 
the others equal to 0. But some of these seven vectors are trivial in the sense 
that their divergence is identically zero, i.e. the conservation equation (13. 4p is 
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satisfied identically. For example, setting in fl3.12p = 1, ai = • • • = 07 = 0, 
i.e. V = y, we obtain the vector 

Ci = 0, C^ = -g{u)uy, C' = g{u)u,, 0^ = 0. 

For this vector Eq. (13.41) is satisfied identically, 

A(C^) + D,{C') + Dy{C^) + D^iC) = -D^{g{u)uy) + Dy{g{u)u,) = 0. 

Let us single out the nontrivial conserved vectors. Since v given by Eq. 
02.81) . the conservation equation fl3.4|) for the vector fl3.12p is written as 

Dt{C') + D,{C') + Dy{C^) + D^iC) = v^F, (3.13) 

where F is given by Eq. (12. ip . 

F = -ut + {f{u)u^)^ + {g{u)uy)y + {h{u)u^);,. 

Then, specifying the expression of Vx from (12. 8p . we write (I3.13P in the form 

DtiC) + DxiC) + Dy{C^) + D,{C^) = [a^yz + a^y + a^z + a,)F. (3.14) 

Eq. (I3.14P shows that we have only four nontrivial conserved vectors. They 
correspond to 01,02,03 and 05, i.e. they are obtained from (13.120 setting by 
turns one of these four parameters to be equal to 1, the others equal to 0. For 
example, the nontrivial conserved vector (13.120 corresponding to 05 is 

= -u, = f{u)ux, = g{u)uy, = h{u)u,. (3.15) 

The conservation equation (13.40 for the vector (13.151) coincides with Eq. (II. ip . 

Thus, the nontrivial conserved vectors are obtained by substituting in 
(13.120 the expression (12. 8p for v with 04 = og = 07 = 0. The resulting vector 

= -{aiyz + a2y + a^z + a^)u, 

= {aiyz + a2y + a^z + a5)f{u)ux (3.16) 
— {aiz + a2)xg{u)uy — {aiy + a3)xh{u)uz, 

= {aiz + a2)g{u){xux + yuy) + {a^z + a^)g{u)Uy, 

= {aiy + a3)h{u){xux + zu-^) + (02?/ + a5)h{u)uz 

is the linear combination with the coefficients 05,01,02,03 of four linearly in- 
dependent vectors, namely the vector (13.150 and the following three vectors: 

= —yzu, = yzf{u)ux — xzg{u)uy — xyh(u)uz, 

= zg{u){xux + yUy), = yh{u){xux + zuz); (3-17) 
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= -yu, = yf{u)u^ - xg{u)uy, 

= g{u){xu^ + yuy), = yh{u)u^] (3.18) 

= —zu, = zf{u)ux — xh{u)Uzi 

= zg{u)Uy, = h{u){xUx + zuz). (3.19) 

The conserved vectors associated with the symmetry X3 from (13. 7p can be 
obtained from the above resuhs merely by the permutation x y, followed 
by the permutations f ■v^ g and ■H- C^. This procedure maps the vector 
(13.121) to the following conserved vector: 

= —UVy, 

= f{u){UyV,+UxVy), (3.20) 

= g{u)uyVy - f{u)u^V:^ - h{u)UzVz, 

= h{u){UyVz + UzVy). 



Accordingly, Eq. (13.141) becomes the following conservation equation for the 
vector (13201): 

A(C^) + D^{C^) + Dy{C^) + D,{C^) = {aixz + a^x + a^z + aQ)F. (3.21) 

It shows that the nontrivial conserved vectors are obtained by substituting in 
(13.201) the expression (12. 8p for v with 03 = 05 = ay = 0. The resulting vector 

C'^ = —{aixz + a2X + a^z + ae)^, 

= {aiz + a2)f{u){xun, + yuy) + {a^z + ae)f{u)ux, (3.22) 

= {aixz + a2X + a^z + aQ)g{u)Uy 

- {aiZ + a2)yf{u)u^ - {aix + a4)yh{u)uz, 

= {aix + a4)h{u){yuy + zuz) + {a2X + aQ)h{u)uz 

is the linear combination with the coefficients of four linearly in- 

dependent vectors, namely the vector (13.151) and the following three vectors: 

= —xzu, = zf{u){xux + yuy), 

= xzg{u)uy — yzf{u)Ux — xyh{u)uz, (3.23) 



C = xh{u){yuy + zUz); 
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= —xu, = f{u){xux + yuy), 

= xg{u)uy - yf{u)u^, = xh{u)uz] (3.24) 

= -zu, = zf{u)u^, (3.25) 
= zg{u)uy — yh(u)uz, = h{u){yuy + zuz). 

Proceeding as above with the symmetry X4 from fl3.7p we obtain, in addi- 
tion to dSjED, flXTTjl - flCT]) and flX^ -f lX^ . the following conserved vectors: 

= -xyu, = yf{u){xu^ + zu^), 

= xg{u){yuy + zuz), (3.26) 

= xyh{u)uz - yzf{u)u^ - xzg{u)uy] 

= —XU, = f{u){xUx + ZUz), 

C'^ = xg{u)uy, = xh{u)uz — zf{u)ux; (3.27) 

C = -yu, = yf{u)ux, (3.28) 
= g{u){yuy + zUz), = yh{u)uz - zg{u)uy. 

Finally, we turn to the time-translational symmetry Xi from (13. 7p . In this 
case W = —Ut- Replacing Ut by the right-hand side of Eq. (11.11) we obtain 
from the first equation (13. 6p : 

C^ = v [Dx{f{u)ux) + Dy{g{u)uy) + Dz{h{u)uz)] . (3.29) 

Now we observe that 

vD.j,{f{u)ux) = [vf{u)ux - VxJ^{u)] , 

where we denote J^{u) = J f{u)du and use the equation Vxx = resulting from 
the representation (12. 8p of v. Transforming likewise two other terms in (13.291) 
we write in the divergent form: 

= Dx [vf{u)ux - VxJ^iu)] + Dy [vg{u)uy - Vyg{u)] + Dz [vh{u)uz - VzT-L{u)] , 

where G{u) = J g{u)du, 'H{u) = J h{u)du. Now we can transfer all terms of 
to the components C^, C^, and obtain = 0. The calculation shows 
that after this transfer we will have = C'^ = = = 0. Hence, Xi does 
not lead to a non-trivial conservation law. 
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Thus, we have proved the following statement. 



Theorem 3.1. The translational symmetries (13. 7p of Eq. (13. ip with arbitrary 
coefficients f{u),g{u), h{u) provide ten linearly independent conserved vectors 

flXT^ . (EHD-dSISl) and (^M-<^M- 

Remark 3.1. Eq. (I2.22p is nonlinearly self-adjoint with the substitution 
(I2.30p containing arbitrary solutions a(t,y), /3(t,y), jityy), cr(t,y) of the ad- 
joint equation (I2.24p to the one-dimensional linear heat equation. Therefore 
the conserved vector constructed by the above procedure for Eq. (I2.22p will 
contain arbitrary solutions of Eq. (I2.24p . 

3.2 Conserved vectors for Equation ( 12.181) 

As mentioned in Section II. 1^ the anisotropic heat equation (II. 2p with an ex- 
ternal source q{u) ^ does not have a conservation form. However, Eq. (II. 2p 
can be rewritten in a conservation form 

Dt{C') + D,{C^) + Dy{C^) =0 

if it is nonlinearly self-adjoint, for example, in the special cases (I2.18P and 
(I2.19P . We will find here the conservation form for Eq. (I2.18P . The calculations 
are similar for Eq. (I2.19p . 

We write Eq. (I2.18P in the form 

ut = f{u)u,j,a: + g{u)uyy + f'{u)ul + g'{u)ul + u'^T{u), w = const., (3.30) 

and have the formal Lagrangian 

C = v[f{u)u^^ + g{u)uyy + f'{u)ul + g'{u)ul + u^J^{u) - ut]. (3.31) 

For this formal Lagrangian Eqs. (II. lip yield (cf. Eqs. (13. 6p ) 

= -Wv, 

C^ = W [f\u)u,v - f{u)v^] + f{u)vD.,{W), (3.32) 

C'' = W [g'{u)uyv - g{u)vy\ + g{u)vDy{W). 

Eq. (I3.30p admits the three-dimensional Lie algebra spanned by the oper- 
ators Xi,X2,X3 from (13. 7p . Let us apply the formula (I3.32p to the symmetry 
X2. In this case W = —u^ and (I3.32p is written (see Eqs. (13. 8p ) 

= vu^, 

= -f'{u)vul + f{u)u^v^ - f{u)vuxx, (3.33) 
= -g\u)vUxUy + g{u)uxVy - g{u)vuxy, 
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where v should be replaced by its expression (l2.2Up . 

V = {Aiy + Bi) cos{ujx) + {A2y + B2) sm{ujx). fl2:20|) 
Let us simplify the vector fl3.33p in the same way as in Section 13.21 We write 

= C^ + D^{uv), 



where 

and replace by 
Hence 



= -uv^, (3.34) 
= C^ + Dt(uv). 



u 



= C^ + vut 

= C \v{f{u)u.,^ + r{u)ul + Dy{(j{u)uy) + u?7{ 

The substitution of the expression fl3.32p for yields: 

= j{u)u^Vx + b?T[u) + vDy[g{u)uy). 

One can verify that the following equation holds: 

vDy[g{u)uy) = Dy [vg{u)uy - Q{u)vy\ . 

It is obtained by introducing the function Q{u) = f g{u)du and noting that 
the equation Vyy = is valid for the representation f l2.20p of v. Hence can 
be reduced to 

= f{u)u^v, + LO^J^{u), (3.35) 

whereas becomes 

C'' = C'' + [vg{u)uy - g{u)vy] . 
The latter equation upon inserting the expression fl3.32p for yields: 

= g{u)uyV^ - Q{u)v,j,y. (3.36) 
Collecting Eqs. (I3.34p - fl3.36p and ignoring the tilde we arrive at the vector 

= -uv^, 

C = f{u)u,v, + co^T{u)v, (3.37) 



= g{u)uyV^ - Q{u)v, 



xy-i 
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where v should be replaced by its expression fl2.20l) . The vector fl3.37p satisfies 
the conservation equation in the following form: 

Dt{C^) + D^{C^) + = V, [{f{u)u,), + {g{u)uy)y + uj^F{u) - ut] . 

The expression fl2.20l) for v contains four arbitrary constants Ai, A2, Bi, B2. 
Accordingly, the vector (13.371) is a linear combination of four linearly indepen- 
dent vectors. Hence, we have demonstrated the following statement. 

Theorem 3.2. The invariance of Eq. (I3.18P with respect to the one-parameter 
group of translations of x with the generator X2 provides the following four 
linearly independent conserved vectors: 

= sm{ujx) u, = — sin(a;x) f{u)ux + OJ cos(a;x) J^{u), 

= — sm{ux) g{u)uy; (3.38) 

= cos{ujx) u, = — cos{ux)f{u)ux — ijj sin(a;x) T{u)^ 

= — cos(u;x) (?('u)'Uy; (3.39) 

= y sin(a;x) u, = —y sin(ci;a;) f{u)ux + uy cos{ijJx) J-'(m), 

= —y sin((X'x) g{u)uy + sin(co'x)^('u); (3.40) 



= y cos{ujx) u, = —y cos{ujx) f{u)ux — coy sin(a;x) J-'(m), 

= —y cos{ijJx) g{u)uy + cos(co'x)^('u). (3-41) 

Remark 3.2. The conserved vectors provided by the generator X3 of the 
group of translations of y can be computed likewise. The generator Xi of the 
time-translations provides only the trivial conserved vector. 

Remark 3.3. Using the conserved vectors (I3.38p - (l3.4ip one can write Eq. 
(12.180 in four different conservation forms. For example, the vector (13.380 
satisfies the conservation equation 

= sin(a;a;) [u^ - {f{u)u,), - {g{u)uy)y - u^J^{u)] . 

Accordingly, Eq. (12.181) can be replaced by the following conservation equation: 

Dt [sm{ijjx) u] — Dx [sm{ijjx) f{u)ux + uj cos(u;x) J^iu)] 

- Dy [sin{ux) g{u)uy] = 0. (3.42) 
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